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ABSTRACT 

Type II SNe can be used as a star formation tracer to probe the metallicity distribution of global 
low-redshift star formation. We present oxygen and iron abundance distributions of type II supernova 
progenitor regions that avoid many previous sources of bias, and can serve as a standard of compari- 
son for properly observationally evaluating how different classes of supernovae depend on progenitor 
metallicity. In contrast to previous supernova host metallicity studies, this sample is homogeneous 
and is drawn from an areal rather than a targeted survey, so supernovae in the lowest-mass galaxies 
are not excluded. We spectroscopically measure the gas-phase oxygen abundance near a represen- 
tative subsample of the hosts of type II supernovae from the first-year Palomar Transient Factory 
(PTF) supernova search. The median metallicity is 12+log(0/H) = 8.65 and the median host galaxy 
stellar mass from fits to SDSS photometry is 10 9 ' 9 Mq. Though iron abundance is more central to 
the evolution of massive stars than oxygen abundance, it cannot be measured directly in extragalactic 
HII regions. Using the relationship between iron and oxygen abundances found for Milky Way disk, 
bulge, and halo stars, we can translate our distribution of type II SN environments as a function of 
oxygen abundance into an estimate of the iron abundance, and find the median [Fe/H]=— 0.60. 
Subject headings: (stars:) supernovae: general; galaxies: abundances; galaxies: dwarf 



1. INTRODUCTION 

The question of how different classes of supernovae 
(SNe) depend on progenitor metallicity has been lim- 
ited by the lack of an unbiased standard of comparison. 
The metallicity distribution of galaxies is insufficient be- 
cause core-collapse supernovae (CCSNe) trace star for- 
mation rather than stellar mass. Existing supernova host 
galaxy metallicity distributions are limited by heteroge- 
neous supernova samples and supernova surveys that do 
not search in the faintest hosts. We will address this by 
creating a metallicity distribution of type II supernova 
progenitor regions that avoids many previous sources of 
bias. 

The observational properties of CCSNe span a broad 
range of spectral types, luminosities, apparent kinetic 
energies, and other properties. Interpreting this diver- 
sity remains a fundamental theoretical and observational 
challenge, particularly as to how differences in the stel- 
lar progenitors of the SNe are related to the explo- 
sions. For example, the relative fractions of hydrogen- 
rich type II SNe and hydrogen-poor type Ib/c SNe can 
be predicted as a function of metallicity b ased on models 
for mass-loss from t he progenit or stars (Eldridg e et al.l 
[20081: iGeorgy et all [20091 |2012| ). Standard mass loss 
mode ls for massive stars are based on line-driven winds 
(e.g. iKudritzki fc Puis! I2000D : the efficiency of these 
winds depends on metallicity because metals, particu- 
larly iron, dominate the line opacities driving the winds 
(jVink fc de Ko tcr 2005) . Furthermore, the explosion en- 
ergy of nor mal type II SNe may dep end on the progenitor 
metallicity (jKasen fc Wooslevll2009D . and certain types of 
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SNe may occur only for low-metallicity progenitor stars 
(lOber et al.l 119831: iHeger fc Wooslevl l2002t lLanger et all 
12007ft . 

Serendipitous observations of SN progenitor stars prior 
to explosion are profoundly useful for sorting out how 
the S N properties depend on the progenitor (e.g. ISmartt] 
120091 and references therein). Unfortunately, observing 
or strongly constraining the properties of the progenitor 
star is only possible for nearby galaxies, leading to very 
small samples. This limits the utility of this technique for 
constraining the properties of subclasses of events that 
are rare and therefore observed mostly at large distances, 
such as extremely optically luminous CCSNe. 

Observational techniques for addressing this question 
without pre-explosion data involve estimating the pro- 
genitor properties from the environments that remain 
behind. These techniques range from estimating pro- 
genitor age a nd mass from the degree o f correlation with 
Ha emission ([Anderson fc Jamesll2008D , to detailed char- 
acterization of resolved stellar populatio ns near super- 
nova explosion sites in n earby galaxies (Badcncs et al.l 
120091 : [Murphy et aLll2011[ ) to extrapolation of stellar pop- 
ulation properties from galaxy photometry and spectro- 
scopic observations (e.g. Kelly fc Kirshneri[20lll ). 

Spectroscopic observations of H II regions near a su- 
pernova can be used to estimate the metallicity using 
strong-line oxygen abundance indicators. There are pre- 
cision and accuracy limitations to strong-line abundance 
estimates, but more rigorous abundance measurements 
using faint auroral lines are too costly (or completely un- 
feasible) for statistical surveys of SN host metallicities. 

The frequency of several classes of CCSN types 
have been found to va ry with host metallicity. 
iPrantzos fc Boissierl ([2003H found that the ratio be- 
tween type Ib/c and type II SNe declines with increas- 
ing host luminosity, and because of the galaxy lumi- 
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nosity/metallicity relationship, they suggested this was 
probably a metalli city effect. This w as confirmed and 
expanded upon by iPrieto et al.l (|2008t ) in a subsequent 
study, which found that hosts of type Ib/c SNe are higher 
metallicity than hosts of type II and type la SNe based 
on spe ctroscopically measured metallicity. Stane k et al.1 
(2006) found that type Ic SNe associated with nearby 
long gamma-ray bursts (GRBs) were in faint, metal- 
poor galaxies, proposing a progenitor m etallicity cut- 
off ab ove which GRBs do not occur, and Modjaz et al. 
(2008) found that type Ic SNe that were associated with 
GRBs were in more metal-poor re gions than those that 
were not. lAnderson et al.l (|2010f ) found that type lb 
and type Ic SNe hosts have marginally higher metallicity 
than hosts of type II and type la SNe. The metallic- 
ity lo£aj_to_tv_pe IcSNe without broad lines was found 
by iModiaz et all (|201lD to be on average higher than 
near type lb SNe, regardless of which strong-line metal- 
licity diagnostic was use d, consis tent with the results o f 
iKellv fc Kirshnerl (|201~TI) . th ough I Anderson et al.l (f2010h 
disagree. IKellv fc Kirshnerl (f201lT found that although 
hosts of type lb and type Ic SNe are higher in metallicity 
than hosts of type II SNe, hosts of broad-lined type Ic 
(Ic-BL) SNe are lower in metallicity. 

Type II SNe which are abnormally optically lumi- 
nous appear to occ ur more often in low-mass, low- 
luminosity galaxies (|Neill et al.l 120111) . Spectroscopic 
abundance measurements of hosts of five such lumi- 
nous CCSNe indic ate this is likely due to metallic- 
ity, as shown by IStoll et al.l (120111) (includi ng data 
from lYoung et all I2010L IKozlowski et al.l l2010h . Sub- 
sequent spectroscopic measurements of the host o f the 
luminous SN 2008am (iChatzopoulos et al.l 120111) and 
the luminous SN 2010av (IPrieto fc Filippenkd l2010l: 
iModjaz et al.ll20Tot iSanders et al.ll201l[ ). and photomet- 
ric limits on th e hosts of the lumino us SNe PSl-lOawh 
and PSl-lOky ()Chomiuk et al.l T201 lh reinforce this con- 
clusion. 

Ma ny studies (e.g. ISanders et al.l 120111: IStoll et all 
120111: iCampisi et all 1201 It iVergani et al.ll2011f ) compare 
host galaxies of SNe or GRBs to the overall galaxy pop- 
ulation, which is a good way to put the host metallicity 
results in context. It is not, however, a secure way to 
evaluate whether metallicity is a key parameter govern- 
ing whether a supernova has a given spectral type or 
luminosity. Supernovae trace star formation rather than 
overall stellar mass, and star formation is not evenly dis- 
tributed among galaxies of a given mass. Rec ent results 
(|Lara-L6pez et al.ll2010l:lMannucci et al.ll2010D show that 
the scatter in the galaxy mass-metallicity relationship is 
reduced by considering star formation as a third param- 
eter, and that star formation rates (SFRs) are higher in 
lower metallicity galaxies at a given mass (but see also 
lYates et al.ll2012f) 

iLevesaue et al I (|2010f ) and Irian etaLl (|2010l ) suggested 
that instead of a metallicity cutoff for GRBs there is 
a separate luminosity-metallicity relationship for GRB 
host galaxies offset to lower metallicity t han the normal 
galaxy mass-metallicity relationship of iTremonti et al.l 
(2004). (These results are consistent with a simple metal- 
licity cutoff with the exception of a single high-mass GRB 
host galaxy that appears to be high metallicity.) This 
question of an offset in the mass-metallicity relationship 
has confused many subsequent discussions of a metallic- 
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iMannucci et al.l (|201lD notes that for galaxies at a 
given mass, lower metall icity galaxies have hig her aver- 
age star formation rates ([Mannucci et al.ll2010l ) and thus 
core-collapse events, which trace star formation, should 
have a mass-metallicity distribution shifted to slightly 
lower metallicity at a given mass . This effect is s imilar 
eno ugh to the o ffset p ropos ed by ILevesaue et al.l (120101 ) 
and lHan et al.l (|20101) that IMannucci et al.l (|201lD con- 



clude that GRB hosts do not differ substantially from 
the typical galaxy population and therefore there is no 
metallicity dependence to GRB hosts. Many subsequent 
studies have claimed to disprove a metallicity connec- 
tion for GRBs or luminous SNe by showing that the 
host galaxy metallicity is consistent with its mass and 
star form ation rate according to the relationship betwee n 
the three (|Lara-L6pez et all 2010: M annucci et al.l l2010). 
While this is evidence against a distinct mass-metallicity 
relationship for the hosts of these SNe, it is not evidence 
against a metallicity dependence. 

These investigations have not yet quantified the ex- 
pected metallicity shift in the star formation rate due 
to this relationship. This shift can be quantified semi- 
analytically by convolving the galaxy mass function with 
the relationship between galaxy mass, metallicity, and 
star formation rate to find the overall di stribu t ion o f 
star formation as a function of metallicity. iNiinol ([201 ll ) 
does so by comparing the observed metallicity distribu- 
tion of GRB hosts to the metallicity distribution of star 
formation, calculated two different ways. First, calcu- 
lated observationally from the stellar masss function, the 
galaxy M-SFR relation, and the galaxy mass-metallicity 
relation. Second, calculated using the relationship be- 
tween galaxy mass, metallicity, and SFR defined by 
IMannucci et al.l (|2010f ) , and assigning metallicities based 
on SFR as well as mass. He finds the difference between 
these two estimates is less than 0.5 dex in oxygen abun- 
dance on the KK04 scale, Regardless which method is 
used, he finds that the GRB host metallicity distribution 
is incompatible with the metallicity distribution of star 
formation unless the GRB fraction depends on metallic- 
ity. (This is not the primary result of that study, which 
examines the fact that galaxies do not have one single 
metallicity, but show an internal spread. Assuming a hy- 
pothetical transient phenomenon which ha s a strict cut- 
off metallicity above which it cannot occur, iNn^lpOllf ) 
shows that such a spread in internal galaxy metallicities 
would serve to widen the observed host metallicity dis- 
tribution of that transient phenomenon.) 

Only by comparing the metallicity distribution of a 
SN variety to the overall metallicity distribution of star 
formation can one rigorously test for a metallicity depen- 
dence. Our metallicity distribution of type II SN hosts 
can serve as a standard of comparison for evaluating how 
different classes of SNe depend on progenitor metallicity. 

A problem for essentially all of these previous studies 
relating metallicity and SN properties is that they draw 
on heterogeneous SN s amples. At pre sent, the Palo- 
mar Transient Factory (|Rau et al.l [20091 PTF) appears 
to supply the closest approximation to an unbiased SN 
sample, in the sense that the SN selection is essentially 
independent of host galaxy properties. This would not 
be true, for example, of the Lick Obser vatory Supernova 
Search (LOSS) survey (|Li et al.l 1201 If) . which explicitly 
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targets larger galaxies. In this paper we will focus on 
the 52 type II SNe found in the first year of PTF op- 
erations. Because some subtypes of CCSNe are known 
to have different distributions in host metallicity, we fo- 
cus on type II SNe. We have measured metallicities for 
the environments of a representative subsample of 34 of 
these type II SNe, and present the resulting metallicity 
distribution. 

This distribution probes the low-redshift star forma- 
tion rate as a function of metallicity in an independent 
way from current methods relying on galaxy population 
statistics. This is an important distribution to charac- 
terize because in order to determine whether a massive 
star outcome has a metallicity dependence, we need to 
examine its frequency relative to the metallicity distribu- 
tion of star formation, not to the metallicity distribution 
of existing stellar mass (the galaxy metallicity distribu- 
tion). 

We place the metallicity distribution of type II SNe en- 
vironments in context by comparing it with previous su- 
pernoya host studies dPrieto et alj[2008t lAnderson et al.l 
mM IKellv fc Kirsh neri 120111). with the SPSS DR7 
MPA/JHU value-added ca talog (iKauffmann et ahl 
20031: iTremonti et al.l 120041 : iBrinchmann et al.l 12001 
Salim et al.l 12007^ an d with estimates f rom galaxy 
population statistics ( Stanek et al.l I2006D . Noting 
that iron is more fundamental to stellar evolutionary 
outcomes than oxygen, we then translate our oxygen 
abundance distribution of type II SNe environments into 
an iron abundance distribution by using the observed 
relationship between oxygen and iron abundance in 
Milky Way bulge, disk, and halo stellar abundances. 

2. SPECTROSCOPIC OBSERVATIONS AND ANALYSIS 

We drew our targets from the fi rst-year of the Pal omar 
Transient Factory (PTF) survey (jArcavi et al.ll2010( ). an 
areal rather than a targeted survey, so supernovae in 
the lowest-mass galaxies are not excluded by selection. 
The sample selection for the survey is not yet published. 
There are 52 type II SNe in the full first-year PTF CCSN 
sample, of which we have measured spectroscopic metal- 
licity determinations for a subsample of 34. This sub- 
sample is characteristic of the overall type II sample, as 
we show in i j4.ll 

We obtained host galaxy spectra of these SNe us- 
ing the Ohio State Multi-Object Spec trograph (OSMOS, 
iMartini et all 120111 : iStoll et aljEoiO: on the 2.4-m Hilt- 
ner telescope, the Wide Field Reimaging CCD Camera 
(WFCCD) on the 2.5-m du Pont telescope, and the dual 
imaging spectrograph (DIS) on the 3.5-m Astrophysi- 
cal Research Consortium telesco pe. We also use twelve 
archival spectra from SPSS PR7 (lAbazaiian et al.ll2009t 
lUomoto et aHHoo! lYork et al.ll200rT1Cunn et al.ll200ffl . 
The properties of the spectroscopic observations are sum- 
marized in Table [Tj Images of nine representative type II 
hosts spanning the observed range of metallicity and 
galaxy mass are shown in Figure [T] We also measured 
spectroscopic metallicities for the hosts of three type lb, 
two type lib, three type Ic, and one type Ic-BL from the 
first-year PTF sample. 

The new observations were made either at the super- 
nova position or at a similar galactocentric radius to min- 
imize any biases from metallicity gradients in the host 
galaxies. We include the angular distance from the host 
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Figure 1. Nine type II SN hosts, spanning the metallicity and 
mass range, arranged so mass increases toward the right and oxygen 
abundance increases toward the top. Each panel is scaled to a 
physical size of 50 kpc (Ho,^ m ,^A = 70,0.3,0.7) and is centered 
on the position of the supernova. All images are from SDSS and 
were taken before the supernovae. 

galaxy center to the supernova site in Table [2] for easy 
comparison to the seeing and to the spectroscopic aper- 
ture, listed in Table[T] We include the projected physical 
distance to facilitate future comparison to studies that 
use galactocentric spectra. For host galaxies with multi- 
ple SPSS spectra or with spectra from multiple sources, 
we fit a metallicity gradient where possible and provide 
the best fit metallicity at the galactocentric radius of the 
supernova progenitor; these are labeled as 'grad'. 

We primarily consider host metallicities dete rmined 
with the N2 diagnostic of iPettini fc Pagell (|2004l ). which 
we directly measure for each of our targets. This di- 
agnostic depends solely on [N II]A6584/HaA6563, and is 
extremely insensitive to reddening, though it has a larger 
intrinsic scatter than other strong-line diagnostics based 
on the physical conditions of the H II regions. There 
are a number of techniques that are used to estimate the 
oxygen abundances of H II regions in star-forming galax- 
ies, and a substantial litera ture discussing their va rious 
merits and drawbacks (e.g. IKewlev fc Ellison! 120081 and 
references therein). A full recap of this is outside the 
scope of this paper, but we discuss the consequences of 
these uncertainties for our study in Section 14.61 

In Table [5] we list the metallicities of the progenitor 
regions of 34 type II SNe, three type lb, two type lib, 
three type Ic, and one type Ic-BL from the PTF first- 
year core-collapse sample. All subsequent analysis is of 
the type II hosts, for which we have good statistics. We 
do not include type lib SNe in the type II sample because 
their spectral similarity to type lb SNe at all but early 
times can lead to some typing issues, and because pre- 
vious studies (e.g. iModiaz et al.ll201lt IKellv fc Kirshneri 
1201 If ) considered them in the stripped-envelope subclass. 
We show the metallicity distribution of the type II hosts 
for several different strong-line metallicity diagnos t ics us - 
ing the empirical conversions of IKewlev fc Ellison! ((20081) 



4 



Stoll et al. 



1 

e0.8 
o 

o 

CO 

to.6 



1 1 1 1 1 1 1 


1 1 | 1 1 1 1 


- PP04N2 




" PP0403N2 


f -J 

ffl if 


; D02 f 




'_ T04 




,- ] ' p 

J Jj 


i Z94 




! 1 j 




1 r 1 i 

r i 


f ■ l J 

r . ; i 


1 1 


r 


LS 


/; ''' ! r 
i fTw it"i i If i i 


I I 1 I I I I 



8.5 9 
12+Iog(0/H) 



Figure 2. Cumulative distribution of local gas-phase oxygen 
abundances for type II SNe. As type II SNe trace young stel- 
lar populations, this traces the metallicity distribution of star for- 
mation at low redshif t. The solid black distribution uses the N2 
diagnostic of Pcttini & Pagcl (2004, PP04) we adopt as our stan- 
dard. Subsequent figures only use these measurements. The other 
distributions show how the result would change for other strong- 
line metallicity di agnostics, based on the empirical conversions of 
IKewlev & 'E llison ( 2001). The dotted orange curve uses the PP04 
Q3N2 diagnostic, the short-dashed gree n curve us es Dcnicol oet al.l 
(2002), the long-dashed blue c urve uses [ Trcmonti ct al. (2004]), the 
dot-dashed purple curve uses IZaritskv et al.l 119941). We do no t 
show conversion s to the scales of Kobulnickv & Kcwlcv l|2004l 'l. 
IMcGaug h (1991), and Kc wlev fc Dopital (120021) . which require ex- 
ternal branch information. 

in Figure [21 In subsequent figures, we ch oose as our scale 
conve ntion only the N2 diagnostic of iPettini fc Pagell 

pool . 

3. CHARACTERIZING THE SPECTROSCOPIC 
SUBSAMPLE 

To investigate whether we have acquired spectra of a 
representative subsample of the hosts of first-year PTF 
CCSNe, we compared the stellar mass and star formation 
rates of the hosts with and without metallicity estimates. 
We used SED models of the SDSS photometry of the 
hosts within the DR8 footprint to estimate the masses, 
SFR, and characteristic stellar ages. Of the 52 type II 
SNe in the full sample, 47 have SDSS photometry. We 
also analyzed the properties of the 19 non-type II PTF 
CCSN hosts that fell in the DR8 footprint, but we will 
restrict our comparisons with our spectroscopic sample 
of type II SNe to these 47 type II SNe to avoid any of the 
currently known selection effects with metallicity linked 
with supernova type (see 

3.1. Extracting fluxes from SDSS imaging 

We obtained ugriz images for the 66 superno va fields 
in the SDSS Data Release 8 (jAihara et aljhoilD . These 
images are fully calibrated in the SDSS natural system, 
which is close to the AB system, and sky-subtracted. 
We combined the most sensitive SDSS bands (gri) for 
each supernova field in order to make a deeper stacked 



image that can be used to find all the galaxies and de- 
fine their photometric apertures. We used these deeper 
stacked images as the re ference image for sourc e de- 
tection using SExtractor ([Bertin fc Arnoutsl I1996D and 
checked by eye the positions around each supernova in 
order to select the most likely host galaxy. We were 
able to assign likely host galaxies for 64/66 SNe. The 
two events without host galaxy detections, PTF09be and 
PTF09gyp, have sources that are > 13 kpc (projected) 
from the positions of the SNe. We note that the host of 
P TF09gyp has a repo rted magnitude of r = 21.75 mag 
in lArcavi et al.l (|2010f) from SDSS photometry, but we 
do not confirm this detection. After selecting the host 
galaxies in the stacked images, we used imedit in IRAF0 
to mask nearby stars, which could contaminate the flux 
measurements, filling the masked regions with the local 
background. Finally, we ran SExtractor on the individual 
ugriz images using the apertures defined from the deeper 
stacked images to obtain total (AUTO) galaxy fluxes. 
W e applied the small (< 0.04 mag) corrections derived 
in iKessler et ail ([20091) to transform the SDSS fluxes to 
the AB system. The resulting coordinates and fluxes of 
the host galaxies are presented in Table [31 and absolute 
magnitudes are presented in Table HI We include 3cr 
upper limits for the hosts of PTF09be and PTF09gyp, 
which were calculated assuming a circular aperture of 
radius r = 5 kpc at the distance of the SN. 

3.2. Galaxy properties 

We used the code for Fitting and Assessment o f Syn- 
thetic Templates (FAST v0.9b. iKriek et alj|2009l) to fit 
these host galaxy spectral energy distributions to esti- 
mate the stellar mass, star formatio n rates (SFRs), and 
charac teristic ages. We chose the iBruzual fc Chariot! 
(2003) libraries, a Salpeter IMF, and Solar {Z = 0.02) 
metallicity to do the fits. We also assumed an exponen- 
tially declining SFR model with r = 1 Gyr for the star 
forming component of the model. The results are pre- 
sented in Table [5j 

Since the star-formation rates derived with FAST can 
be dependent on the assumed SFR model, we obtained 
an independent SFR estimate using the M-ba nd luminosi- 
ties of the galaxies. We used the results of iSalim et al.l 
(2007) for - 50000 SDSS galaxies with GALEX pho- 
tometry to derive a relation between the absolute it- 
band magnitudes, corrected for intrinsic attenuation, 
and their derived SFRs. We obtain a linear relation 
fit between M u and SFR ( in M Q /yr) of the form 
log(SFR) = -0.36 xM u -6.73 (for a Salpeter IMF), valid 
for 2 > log(SFR) > -2 and log(SFR/M*) > -10.5, with 
an rms scatter of 0.23 dex. We applied this relation to ob- 
tain SFRs for the supernova h ost galaxies, u sing t he Low- 
Resolution Template code of lAssef et al.l (|2008| ) to de- 
rive -ftT-corrected it-band absolute magnitudes corrected 
for Galactic extinction. After further correcting these 
magnitudes by intrinsic attenuation using the values ob- 
tained with FAST and the Calzetti redden ing law, we 
applie d the linear relation derived from the ISalim et al.l 
(|2007h data to estimate SFRs. These independent val- 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Figure 3. The two methods of estimatin g SFR are c onsistent. 
SFR estimates based on u-band photometry (Salim ct al. 2007) arc 
plotted against SFR estimates based on FAST for the full sample 
of PTF type II SN hosts (open red points) and the metallicity 
subsample (solid black points). The blue dashed line represents a 
1:1 correspondence. 

ues are presented in Table [5] The agreement with the 
SFRs derived by FAST is fairly good in general, with 
a Kolmogorov-Smirnov (K-S) test probability of 81% of 
the results of the two different methods being drawn from 
the same underlying distribution. The two SFR estima- 
tion methods are directly compared in Figure |31 The 
SFR calculated from the u-band luminosity with aper- 
ture corrections is more consistent with the method the 
MPA/JHU value-added catalog uses to determine star 
formation rates than the FAST template fitting. 

4. DISCUSSION 

This is the largest homogeneous sample of supernova 
host galaxy spectra metallicity measurements yet. Nev- 
ertheless, limited observing time made following up the 
entire PTF type II sample impractical, so we first con- 
sider whether our spectroscopic host sample is represen- 
tative of the full sample, finding that it is. We then 
place the SN hosts in context by comparing their prop- 
erties to those of galaxies in the MPA/JHU value-added 
catalog. While they are well-matched in galaxy mass and 
metallicity, the type II hosts appear to be biased toward 
higher star formation rates than the galaxies in the cata- 
log. The metallicity distribution of these type II hosts is 
remarkably similar to that found by previous studies of 
hosts of type II SNe, despite their more inhomogeneous 
selection. Their metallicity distribution is also consis- 
tent with a distribution of star formation calculated from 
galaxy population statistics. We discuss how our study 
avoids some selection effects due to supernova type and 
host galaxy type that might influence the metallicity dis- 
tribution. A key future use of the metallicity distribution 
of type II SNe we find here will be to evaluate possible 
metallicity dependence of other subclasses of CCSNe. We 
discuss the advantages and disadvantages of the metallic- 
ity diagnostic we choose for this study. Finally, we define 



a relationship between oxygen and iron abundances, and 
convert our observed oxygen abundance distribution into 
an assumed iron abundance distribution, as iron is more 
important to the evolution of massive stars than oxygen. 

4.1. How representative is the metallicity sample? 

To investigate how representative the subsample for 
which we have spectra and metallicities is of the entire 
sample of PTF type II SN hosts we compared the distri- 
butions of the two samples in galaxy mass, characteristic 
stellar age, and star formation rate, using the 32 (47) 
hosts with (without) measured metallicities that also lie 
in the SDSS DR8 imaging footprint. Figures 0] and [5] 
show that the two sub-samples have essentially identi- 
cal distributions in host mass (K-S probability 99.9%), 
age (57.8%), and SFR (64-74%, depending on the SFR 
estimation method). 

We also investigated the effects of redshift on the com- 
pleteness of the sample by dividing it into lower and 
higher redshift subsamples and comparing the proper- 
ties of the two. In Figure [5] we show the metallicity 
distribution of these two subsamples. A K-S test indi- 
cates that the two have a 19% probability of being drawn 
from the same metallicity distribution, consistent at ap- 
proximately la. The hosts in the two redshift bins have 
essentially identical distributions in host mass (K-S prob- 
ability 63%), age (91%), and SFR (91-99%, depending 
on the estimation method). 

4.2. The type II metallicity sample in context 

We next compare our sample to galaxy properties 
in the D R7 SPSS MPA/JHU value-added catalog for 
context (iKauffmann et all I2003L ITremonti et al.1 120041 : 
iBrinchmann et alJl2004t iSalim et al.ll2007| ). We compare 
to the subset of the DR7 objects which have redshifts 
within the range of our sample, successful estimates of 
the stellar mass and star formation rate, a 12+log(0/H) 
metallicity estimate, and an [N II]A6584/HaA6563 flux 
ratio within the valid range for the PP04 N2 metallicity 
diagnostic. This last requirement has almost no effect, 
reducing the sample by only 0.7%. The strictest condi- 
tion b y far is the requirement of a valid ITremonti et al.l 
(2004) metallicity estimate. Note that the galaxy prop- 
erties in the MPA/JHU value-added catalog are a func- 
tion of redshift because is not a volume-limited sample. 
The selection that defines the catalog is reflected in the 
properties of its constituents. 

As shown in Figure [3 the SN hosts appear to trace the 
MPA/JHU sample well. They do not appear to be biased 
toward lower metallicities at a given mas s (left panel) as a 
simpli stic interpretation of the results of Man nucci et al.l 
(2010) might suggest. It has been suggested that core- 
collapse SNe and GRBs may be less frequently observed 
in higher m etallicity environments due to higher ex- 
tinction (e.g iMaiolino et al.ll2002[: iMannucci et alj|2003t 
iCresci et alll2007t iCampisi et al.N2011| ). but this sample 
does not show evidence for such a bias. The SN hosts 
do, however, appear to be biased toward higher star for- 
mation rate^l (right panel) than the MPA/JHU galaxy 
sample, as would be expected if they do indeed trace star 
formation. They appear to well-sample the distribution 

6 Here we use the SFR calculated from the u-band luminosity 
with aperture corrections. 
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Figure 4. The distributions of the full PTF type II sample (dotted red) and those for which we obtained metallicities (solid black). 
Section 13.21 describes how the host properties were estimated. The subsample for which we have spectroscopic metallicity measurements 
is quite representative of the full sample, with K-S test probabilities of 99.9% (mass) and 57.8% (age) that they are drawn from the same 
distribution. 
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Figure 5. The subsample of type II hosts for which we have measured metallicities (soli d black) is representative of the full sample of 
PTF type II SN hosts (dotted red) in SFR based on u-band photometry l lSalim et al . 2007, left) and FAST (right). 



of the MPA/JHU sample in galaxy mass and star forma- 
tion rate, as shown in Figure [5J 

4.3. Comparing with other SNe host samples 

Work on the m e tallici ty distr ibution of supernova hosts 
by iPrieto et all p00§) and iKellv fc Kirshnerl (|20T1 
looked at overall galaxy metallicity with serendipitous 
SPSS spectra, without isolating the SN site. Studies by 
lAnderson et al.l (|2010f ) tried to measure abundances at 
the SN site or at a similar galactocentric radius, as we 
have done in this study. These previous studies had uni- 
form spectroscopy but the source SN samples were het- 
erogeneous, including SNe discovered in a wide variety 



of ways with very different selection effects. Our source 
SN sample is homogeneous, from a single survey with 
uniform selection. The source survey is areal rather than 
galaxy-targeted, which enables the detection of events in 
the lowest-mass galaxies, removing or at least mitigat- 
ing a possible bias towards high metallicity environments 
which we expect exists in prior supernova host samples. 
We do not attempt to define a volume-limited supernova 
sample here. Rather, we point out that this existing sam- 
ple represents a substantial step forward for this purpose 
fr om previous simila r samp les. 

IKellv fc Kirshner (1201 ll) used abundances follow- 
ing iTremonti et al.l (|2004f > and the 03N2 method of 
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Figure 6. The metallicity distribution of the type II hosts separated into two bins in redshift of seventeen hosts each (left). The green 
solid (red dotted) line is the lower (higher) redshift half of the sample, 0.019 < 2 < 0.035 (0.037 < z < 0.11). On the right, the distribution 
is expanded and metallicities are plotted against redshift. Here the lower redshift half of the sample is green solid points, and the higher 
redshift is red open points. A K-S test shows the two are consistent with being drawn from the same distribution in metallicity at the 19% 
level. There do not appear to be any major selection effects with redshift. 
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Figure 7. Measured host metallicities (black points) as a function of galaxy mass (left) or SFR (right), overlaid on the distribution in the 
MPA/JHU value-added catalog in the same redshift range (blue contours). Our hosts are slightly offset to higher star formation rates, as 
we would expect if type II SNe trace star formation. 

IPettini fc Paeel {2004) turns over at low metallicities, 
and does not match the reverse conversion. To ensure 
monotonicity, we define an ad hoc co nversion by fitting 
to the reverse conversion defined by iKewlev fc Ellison] 
(2008') (see Appendix). This conversion is a fit to a 
fit rather than a fit to data, but it achieves two pur- 
poses: avoiding a turnover at low metallicity and match- 
ing the inverse conversion from PP04 N2 to T04. We use 
this conversion everywhere in this paper when converting 
from the scale of the T04 diagnostic to the scale of the 
PP04 N2 method. Approximately 9% of the lPrieto et al.l 



IPettini fc Page! (l200l. In th is pap er we use the 
N2 method of IPettini fc Pagell (|2004l ). so we must 
convert these to a common scale. The valid range 
of the conversion defined by IKewlev fc Elhsonl (2008) 
from the scale of the 03N2 method to the scale of 
the N2 method does no t span the abundance s here, 
so we convert from the iTremonti et all (|2004| ) abun- 
dances instead. iPrieto et al.l (I2008D al s o uses abundances 
from the method o f ITremonti et aT] (120041). The con- 
version defined by IKewlev fc Ellison! (120081) from the 
scale of ITremonti et all (|2004[ ) to the N2 method of 
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Figure 8. Host mass and SFR for our type II SN hosts with Z 
measurements (large black points), all first-year PTF type II SNe 
hosts in the SDSS photometric footprint (small red points), over- 
laid on the distribution in the MPA/JHU catalog (blue contours). 



(2008) sample have T04 metallicitie s above or below 
the va lid conversion range defined by IKewlev fc Ellison] 
(2008). Because our conversion is well-behaved at low 
metallicity, we convert the entire sample instead of ex- 
cluding that 9%. The decision to include these hosts 
with the lowest and highest metallicity is not highly con- 
sequential; the distribution of the full sample with our 
conversion is consistent (at a K-S probability of 66%) 
with the remaining 91% converted using the method of 
IKewlev fc Ellisor] poM ). 

We compare our type II host metall i city distribution 
with t hose f ound bylPrieto et al.l (120081 ). [Anderson et all 
([2010]) and IKellv fc Kirshnerl (120111) in Figure ! The 
K-S test probabilities that our type II sample could be 
selected from the same underlying distribution as those 
in the earlier studies are 87%, 45%, and 11%, respec- 
tively. The agreement of our SN metallicity distribution 
with the results of these prior studies is striking, given 
the different sample selection. 

4.4. Comparing with SFR metallicity distributions from 
galaxy population statistics 

One semi-observational way of determining the global 
distribution of star formation as a function of metallicity 
is to combine the observed galaxy mass function, the ob- 
served mean star formation rate as a function of galaxy 
mass, and the observed galaxy mass-metallicity relation- 
ship. In Figure [TU] we plot our distribution against one 
such determination based on galaxy population statistics 
(|Stanek et al.| [2006|) Each of the three observed relation- 
ships that are inputs to this alternate method of find- 
ing the metallicity distribution of star formation has its 
own set of selection effects. Potential biases or redshift- 
dependent effects in the samples used to define the re- 
lation could offset the distribution in metallicity. The 
width of the distribution may be misrepresented by us- 
ing a mean relationship to translate from one prop- 



Figure 9. Distribution of local gas-phase oxygen abundance of 
our sample of type II SNe (black solid line) compared with existing 
SN host metallicity samples. The dotted orange curve s hows the 
heterogeneous type II subsample of Anderson et al. ( 2010). The K- 
S probability of this sample being drawn from the same underlying 
distribution as our PTF type II sample is 11%. Similarly, the green 
long- dash curve shows th e distribution of type II, IIP, and Iln SNe 
from IPrieto et aT] 120081 1 (K-S probability 87% ) and the purple 
dash-d ot curve shows the type II sample from Kelly & Kirshnerl 
(2011) (K-S probability 45%). In the latter two cases we converted 
their T04 scale metallicity to the PP04N2 scale using the conversion 
given in the Appendix. 

erty i nto another, such as t he mass-metallicity relation- 
ship (jTremonti et al.l [2004) or the three-way relation- 
ship between mass, metallicity, and star form ation rate 
(jLara-Lopez et ail 12010; M annucci et al.l l2010'). because 
the scatter around that relationship may not be carried 
through to the final distribution. 

Core-collapse SNe, as the deaths of massive, young 
stars, are a relatively good tracer of SFR. Using the 
metallicity distribution of a uniform sample of type II 
SN sites to approximate the metallicity distribution of 
star formation, as we have done, should have almost 
completely independent selection effects from methods 
relying on galaxy population statistics. 

4.5. Selection effects 

One of the primary potential sources of incompleteness 
in using type II SNe environments to trace the metallic- 
ity distribution of global star formation will of course 
be the selection of the sample of type II SNe. An ideal 
sample for this purpose would be a complete, volume- 
limited sample, monitoring a fixed region of sky for a 
fixed period of time, and then eliminating events outside 
the complete sample. Up until very recently, most su- 
pernova surveys have monitored large, luminous galaxies 
rather than regions of the sky, a methodology which has 
the potential to miss any SNe in the very lowest end of 
the galaxy luminosity function. The Palomar Transient 
Factory survey is areal, which removes the potential bias 
against extremely low-mass host galaxies of targeted sur- 
veys. Because the survey selection is not yet published, 
however, we are unable to correct for any biases in our 
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Figure 10. The distribution of star formation as a function of 
metallicity from galaxy population statistics from IStanek et al.l 
(2006) (green dotted line) to our distribution of site metallicities 
of type II CCSNe (black solid line). 

source sample to get a distribution for a volume-limited 
sample. With the current rapid growth in depth and 
breadth of supernova surveys, there is great potential for 
this method of determining the metallicity distribution 
of star formation. 

Another potential source of bias for this method is any 
dependence of the likelihood of a massive star resulting 
in a type II SN on metallicity. Type II-P SNe make 
up around 70% of all type II SNe in the LOS S survey, 
which focuses on relatively luminous galaxies ([Li et al.l 
1201 If ). The frequency of type II SNe has not yet been 
found to depend on metallicity, unlike type lb and Ic SNe. 
We examine only the distribution of type II hosts, for 
which we have good statistics. The type II distribution is 
consistent with the distribution of the nine type Ib/Ic/IIb 
hosts we have measured with a K-S probability of 23%, 
and the distribution of the hosts of all CCSNe we have 
measured (including all the type II hosts) is consistent 
with the type II distribution with a K-S probability of 
99%, shown in Figurc[TT](left). Completeness corrections 
to translate between type II hosts and all CCSN hosts 
may depend slightly on metallicity. 

Several types of CCSNe are known to vary in fre- 
quency with metallicity, as discussed in fjl] Our very 
small sample of spectroscopic metallicity measurements 
of hosts of type Ilb/Ib/Ic SNe is statistically consis- 
tent with these previous results (IModiaz et al.l 120111 : 
lAnderson etldl l201fl Kelly k Kirshnerl 1201 1[) .but not 
at very high significance level. Using the host galaxy 
properties we fit from SDSS photometry gives us slightly 
better statistics, shown in Figure [TT] (right). The relative 
distributions of type II, lib, lb, Ic, and Ic-BL in photo- 
metrically calcu lated host mass ar e cons istent with the 
recent results of iKellv k Kirshnerl (|2011[ ). The overrep- 
resentation of types lib and Ic-BL in low-mass h osts is 
consistent with the results of lArcavi et al.l (|2010l) based 
on host galaxy luminosities. 



4.6. Strong-line metallicity diagnostics 

There is a substantial literature on the merits and dis- 
advantages of each commonly-used method of determin- 
ing metallicities based on fluxes of strong emis sion lines 
(e.g. iBerg et al.l [20TT1: iKewlev k Ellison! 120081 . and ref- 
erences therein). These methods all rely on simplifying 
assumptions about the H II regions being examined: uni- 
formity of electron density, cooling dominated by oxygen 
(implying that other cooling species have abundances 
that vary in lockstep w ith oxygen), and ionization- 
bounded H II regions ?e.g. lPageret^fl979l ). The meth- 
ods can be classified into rough categories: direct meth- 
ods, which rely on estimates of the electron temperature 
and require measurements of faint auroral lines such as 
[O III] A4363 A, empirical (e.g.lPettini k Pagell2004D . the- 
oretical (e.g. iKobulnickv k Kewlevl 120041 ) and acombi- 
nation of empirical and theoretical (e.g. iDenicolo et al.l 
2002). All these methods are based on high S/N mea- 
surements of the line ratios of different combinations of 
emission lines from ions present in the optical region of 
the spectrum (~ 3700 - 6800 A: [O II], [O III], [N II], 
[SII], Ha, and H/3). 

The slope and intercept of the galaxy mass- 
metallicity relationship is different for each diagnostic 
(jKewlev k Ellison! 120081 ). which is a relatively straight- 
forward way to demonstrate that they are not all directly 
measuring some platonic ideal of a fundamental oxygen 
abundance measurement. (This problem is independent 
of the separate question of the exact value of the solar 
oxygen abundance, which also affects how measured gas- 
phase abundances map to stellar abundances.) Instead, 
each technique measures something that correlates well 
with oxygen abundance, but does not directly map to it. 
The simplifying assumptions that allow us to use each 
strong-line indicator to estimate oxygen abundance are 
not perfect for all H II regions. Rigorously selecting the 
best diagnostic for a given situation requires better data 
than are achievable for distant and faint targets, and 
doing a case-by-case selection of strong-line method on 
insufficient data would introduce its own biases. 

The primary advantage of the strong-line techniques is 
that they are possible with fewer photons, and are there- 
fore feasible to perform on large samples for good pop- 
ulation statistics. Given the scale differences between 
methods, however, it is crucially important to ensure 
that all metallicities one is comparing are on the same 
scale. Where possible, we do this by natively determining 
the metallicities in a common scale. Where impossible, 
we convert a metallicity determ ined on another scale us - 
ing the empirical conversions of lKewlev k Ellison! ([2008) 
(but see Appendix [AJ. We emphasize that any meaning- 
ful conclusion should not be affected by a change in the 
metallicity calibration. 

The pri mary scale we choose for this study is the N2 di- 
agnostic of lPettini k Pagel ([2004D . which depends solely 
on [N II]A6584/HaA6563. A disadvantage of the method 
is that it shows a larger dispersion from the average fit 
than for most other strong line diagnostics. There are a 
number of advantages of this diagnostic, however. The 
ratio is very insensitive to reddening due to the close 
wavelength proximity of the lines. It is monotonic with 
oxygen abundance. It depends on lines with relatively 
high fluxes in star forming environments, which means 
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Figure 11. The metallicity distribution (left) of type II SN hosts (solid black), all other CCSN hosts (dashed orange) and all CCSN 
hosts in Table [2] (dotted green). With this small sample of other subtypes, the type II distribution is statistically consistent with the other 
CCSNe (K-S probability 23%) and with all CCSNe (K-S probability 99%). Increasing the sample size, the distribution in photometrically 
calculated host galaxy mass (right) of type II SN hosts (black solid), Ic-BL (pink dotted), lib (green short-dashed), lb (blue long-dashed), 
and Ic (red dot-dashed). 

that good metallicity estimates can be ach i eved at rela- 
tively low observational expense. lYin et abl ()2007f ) find it 
is more c onsistent with T ? metho d s than the 03N2 dia ; 
nostic of iPettini fe Pagell (|2004fl . iBresolin et ail (f20CF 
compare a variety of strong-line abundance estimators 
and find that PP04 N2 is the closest match in both slope 
and normalization to the oxygen abundance gradient in 
NGC 300 measured with stellar metallicity (blue super- 
giants). 



119971: 1 Johnson et all 120071 : lEpstein et al.1 120101 ) (but see 
e.g. llsraelian et al.lll998[ K 

To estimate the conversion, we compared [O/Fe] to 
[Fc/H] over a wide range in metallicities using stellar 
abundance mea surements from the Milky Way bulge , 
disk, and halo (iFulbright et al.l [20071 iRich et all 120071: 

10 



4.7. Iron abundances 

Iron is more fundamentally important than oxygen for 
the late-stage evolution of massive stars, because iron 
provides much of the opacity for radiation-driven stel- 
lar winds (|Pauldrach et al.|[l985 IVink fc de Koten[20q5l 
e.g.). Unfortunately, gas-phase iron abundances are dif- 
ficult to measure, and the fraction of iron depleted onto 
grains is highly variable. Even within our o wn galaxy, 
measuring iron abundances is challenging (e.g. iRodriguezl 
120021 : Pensen fe Snowl[2007l : lOkada et al.ll2008f ) . We mea- 
sure gas-phase oxygen abundances instead as a proxy for 
metallicity because it is observationally feasible. 

Gas-phase oxygen abundance estimated with PP04 
N2 is well correlate d with stellar oxygen abundances 
(jBresolin e t al. 2009). Our gas-phase oxygen abundances 
are thus a good proxy for the stellar oxygen abundances 
of the progenitors of these CCSNe. We can use our oxy- 
gen abundances to estimate the more physically impor- 
tant (for massive star evolution) iron abundances. This 
mapping is not linear because iron abundance varies 
more steeply than oxygen abundance in stars. At low 
metallicity, a-elements like oxygen are enhanced relative 
to iron compared to the solar mixture. In the galactic 
disk and halo, at [Fe/H] > —1, [O/Fe] is approximately 
inversely proportional to [Fe/H], while below [Fe/H] = 
— 1, [O/Fe] may flatten out at a constant (and lower) 
relative iron abundance (e.g iTinslevi 119791 iMcWilliaml 



Rich fc Origliall2005t iLecureur et al.l 12007b iReddv et al.1 
2003L 120061 : iBensbv et al.ll2004UChen et al.ll2003f l Noting 
that [O/H] = [O/Fe] + [Fe/H], we fit the relationship 
between [O/H] and [Fe/H] with an unweighted linear fit, 
as seen in Figure 1121 and find 



[Fe/H] = Cl +c 2 ([O/H]), 



(1) 



where ci = -0.34 ±0.01 and c 2 = 1.25 ±0.05. Although 
the eye is drawn to a steeper trend at higher metallicity 
than the formal fit shown, this is an illusion based on 
relatively few points; the fit is driven to be flatter by a 
dense concentration of points with —0.2 < [O/H] < 0.2 
and —0.5 < [Fe/H] < 0. The relationship does not 
differ substantially between bulge stars (shown in red) 
and halo and disk stars (shown in blue). If [O/Fe] flat- 
tens out below [Fe/H] = —1, the linear relationship we 
choose to fit may not be ideal. We would expect to find 
a slightly steeper relationship were we to exclude points 
below [Fe/H] = — 1, which would mean that type II SN 
progenitors in low oxygen abundance host regions have 
even lower iron abundance than we find here. Because 
there is finite scatter in the measured relationship, how- 
ever, imposing a strict cut at [Fe/H] = —1 actually drives 
the fit to be slightly flatter by biasing the points with 
lowest [O/H] to higher [Fe/H]. 

Applying the fit to a direct conversion between 
12±log(0/H) and [Fc/H] requires assuming a solar oxy- 
gen abundance because stellar abundances are measured 
relative to solar, while gas phase abundances are abso- 
lute (at least to the accuracy of the assumptions for a 
given strong line method). There is currently some dis- 
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Figure 12. To translate oxygen abundances into iron abundances, we fit a linear relation (black) to the iron and oxygen abundances 
of Milky Way bulge, disk, and halo stars. At low metallicity, ct-elements like oxygen are enhanced relative to iron compared to the solar 
mixture. O n the left is the well-known relationship between [O/Fe] and [Fe/H], The red open points are iron and oxygen abundances of 
bulge stars (Fulbright ct al. 2007; Lccurcur ct al. 2007; Rich & Origlia 2005; Rich ct al. 2007). The blue solid points are halo and disk 
stars l IBensbv et aljl2004l:IChen et al.H2003: IReddv et alJI200a. 12009) On the right we express this in terms of [O/H] and [Fe/H] and fit the 
relation. Although the eye is drawn to a steeper trend at higher metallicity than the formal fit, this is an illusion based on relatively few 
points; the fit is driven to be flatter by a dense concentration of points with —0.2 < [O/H] < 0.2 and —0.5 < [Fe/H] < 0. Equation [2] can 
be used to conservatively convert a measured gas-phase oxygen abundance to iron abundance for a given solar oxygen abundance (modulo 
the uncertainties of equating gas-phase strong-line oxygen abundance indicators to stellar oxygen abundances). 



pute over the solar abundance because results from at- 
mospheric and interiors methods differ. For a given solar 
oxygen abundance Oq, 

[Fe/H] = a - c 2 O + ca(12 + log(0/H)). (2) 

Using a solar oxygen abunda nce of O = 8.86 
pelahave fc PinsonneauTtl 120061) . the conversion is 
[Fe/H] = -1 1.4 + 1.25(12 + togf O/H)), while using 
O = 8.69 (|Asplund et al.l I2009D . the conversion is 
[Fe/H] = -11.2 + 1.25(12 + log(0/H)). 

One more assumption must be made; we must convert 
strong-line oxygen abundance to stellar oxygen abun- 
dance. Which strong line indicator maps most precisely 
to stellar oxygen abundance across a wide variety H II 
region conditions in a wide variety of galaxy types is not 
well known. For the purpose of this paper, however, we 
assume that PP04 N2 maps precisely to stellar oxygen 
abun dances (as motivated by the results of lBresolin et al.l 
2009), and make no correction to either slope or zero 
point. 

We apply this fit to transform our oxygen abun- 
dance distribution of type II progenitors into an iron 
abundance distribution, shown in Figure 1131 The me- 
dian value of [Fe/H] is —0-60 u sing the solar value of 
iDelahave fc Pinsonneau lt (2006). If another solar oxy- 
gen abundance O© is assumed, the calculated iron value 
shifts by 

c 2 (8.86-O ), (3) 

so using O = 8.69 (jAsplund et al.ll2009l) . for example, 
the median value of [Fe/H] is —0.39. 

A striking outcome of this translation is that all of the 
type II SN progenitors in this sample appear to have sub- 
solar iron abundances. Although this is notable, it is not 
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Figure 13. The estimated [Fe/H] distribution of the type II SN 
sites. The blue solid line i s the distribution assuming the solar 
oxygen abundance is 8.86 (Dclahavc & Pinsonneault 2006), and 
the red dashed lin e is the distribution assuming the solar oxygen 
abundance is 8.69 HAsplund et al.ir2009h . 

entirely surprising; the sun is more enhanced in iron than 
most Milky Way stars at its oxygen abundance, and all 
of the host galaxies in this sample are smaller than the 
Milky Way. Were we to assume that [O/Fe] flattens out 
below [Fe/H] = — 1 instead of remaining linear, the slope 
we fit would be steeper and the iron abundances we find 
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here would be even lower. 

The type II progenitors all have iron abundances 
greater than [Fe/H] = —1.5, putting them squarely in the 
regime where winds are primarily driven by iron opacity. 
For the most metal-poor stars (Z/Z Q < 1CP 3 ), non-iron 
elements such as carbon dominate the radiative driving 
(|Vink fc de Koterl 12005). but in the metallicity range of 
these type II SN progenitors, iron abundance should still 
be the dominant factor which determines wind strength 
and mass loss. 

5. CONCLUSIONS 

The primary result of this paper is a new progenitor 
region metallicity distribution for a uniform sample of 
type II SNe. Understanding this distribution is impor- 
tant for understanding any possible metallicity depen- 
dences of different types of events associated with mas- 
sive stars, and it can serve as a probe of the metallicity 
distribution of star formation. 

The host galaxies of our type II sample appear to well 
trace galaxies from the MPA/JHU value-added catalog 
in mass and metallicity, showing a slight bias towards 
higher star formation rates. 

We find a similarity between the existing host metal- 
licity distributions for heterogeneous type II supernova 
samples and the metallicity distribution we derive. Be- 
cause the existing host metallicity distributions are based 
on supernova samples that are drawn predominantly 
from galaxy-targeted supernova searches, one might 
naively expect these previous distributions might be bi- 
ased towards higher mass and therefore higher metallicity 
galaxies. We do not find such a trend. 

Comparing to the metallicity distribution of star for- 
mation rather than to the metallicity distribution of 
galaxies as a function of mass is the correct way to evalu- 
ate a possible metallicity dependence of a transient popu- 
lation associated with young stars. We point out that us- 
ing CCSNe to trace star formation leads to an almost en- 
tirely independent way of probing the metallicity distri- 
bution of star formation from methods involving galaxy 
population statistics, and we compare the metallicity dis- 
tribution we derive to one of these. 

Finally, we present our host metallicity distribution in 
terms of iron abundance, by converting our oxygen abun- 
dance distribution to an iron abundance distribution us- 
ing the a/Fe relationship observed in Milky Way bulge, 
disk, and halo stars, noting that iron is more important 
than oxygen for the late-stage evolution of massive stars. 



We find that -1.2 < [Fe/H] < for these type II SN 
progenitors. Though all have sub-solar iron abundance, 
none are metal-poor enough that elements other than 
iron will dominate the wind-driving opacity of the pro- 
genitor star. 

Future improvements to this estimate of the metallicity 
distribution of type II SNe can be made by performing 
completeness corrections for any selection or followup bi- 
ases in the source survey. If the peak luminosity of type II 
SNe is found to depend on host galaxy metallicity, there 
may also be Malmquist-like biases to correct. 
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APPENDIX 
METALLICITY CONVERSIONS 

It is well known that the various strong-line oxygen abun dance estimators have d ifferent scales a nd zero-points. 
(For e xcellent pictorial representations of this, see Figure 2 of IKewlev fc Ellison! (l2008h or Fi gure 12 of iBresolin et al.l 
(2009).) Because of these differences in scale a nd zero-point, it is critic ally necessary to put different estimates on 
a common scale before comparing abundances. IKewlev fc Ellisonl (|2008f ) determined empirical conversions between 
many of the commonly used scales by fitting the trend defined by performing a given two diagnostics on a large sample 
of high S/N galaxy spectra from SDSS. While these conversions have been very useful to the community, there can 
be problems using them at very low metallicity. The forward and reverse conversions between two methods are not 
always consistent, as can be seen in Figure Q31 The problem appears to be a consequence of the interaction of the vast 
statistical weight of the abundant high-metallicity galaxies and the third-order polynomials (which are not precisely 
invertible) in which the conversions are expressed. The high-metallicity end is tightly pinned, allowing low-metallicity 
end of the third-order polynomial to shift between the forward and reverse conversions. 
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Fig ure 14. Forward and rever se conversion over the entire valid conversion range s between the metallicity scale of Zaritskv ct al. (1994) 
and IMcGaugbl QME|) (left) and ITremonti et al.l J2004T) and IPettini fc Paeell J200W) (right). To Z94 and T04 is the red dotted line, while 
from Z94 and T04 is the black solid line. The blue d ashed line would represe nt a conversion between two exactly equivalent metallicity 
scales. The scales are the same as those in Figure 3 in Kcwlcy & Ellison (2008). 

For the purposes of this paper, we need to avoid the double- value of the T04 to PP04N2 conversion at low metallicities. 
We do so by defining an ad hoc conversion from T04 to PP04N2 , as shown in FigurelTSl T his conversion is a third-order 
polynomial fit to the conversion from T04 to PP04N2 given in iKewlev fc Ellison! 1)20081 ). This fit is defined by 

y = a + bx + ex 2 + dx 3 , (Al) 

where y = 12+log(O/H)p P0 4N2, x = 12+log(0/H) T o4, a = 178.248, b = -59.2077, c = 6.80078, and d = -0.257326. 
Note that this is a fit to a fit rather than a fit to data. We define it solely for the purposes of avoiding a turnover at 
low metallicity and matching the inverse conversion from PP04N2 to T04. 



REFERENCES 



Abazajian, K. N, et al. 2009, ApJS, 182, 543 
Aihara, H., et al. 2011, ApJS, 193, 29 

Anderson, J. P., Covarrubias, R. A., James, P. A., Hamuy, M., & Habergham, S. M. 2010, MNRAS, 407, 2660 
Anderson, J. P., & James, P. A. 2008, MNRAS, 390, 1527 
Arcavi, I., et al. 2010, ApJ, 721, 777 

Asplund, M., Grevesse, N, Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Assef, R. J., et al. 2008, ApJ, 676, 286 

Badenes, C, Harris, J., Zaritsky, D., & Prieto, J. L. 2009, ApJ, 700, 727 
Bensby, T., Feltzing, S., & Lundstrom, I. 2004, A&A, 415, 155 
Berg, D. A., Skillman, E. D., & Marble, A. R. 2011, ApJ, 738, 2 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 

Bresolin, F., Gieren, W., Kudritzki, R., Pietrzyflski, G., Urbaneja, M. A., & Carraro, G. 2009, ApJ, 700, 309 

Brinchmann, J., Chariot, S., White, S. D. M., Tremonti, C, Kauffmann, G., Heckman, T., & Brinkmann, J. 2004, MNRAS, 351, 1151 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Campisi, M. A., Tapparello, C, Salvaterra, R., Mannucci, F., & Colpi, M. 2011, MNRAS, 417, 1013 
Chatzopoulos, E., et al. 2011, ApJ, 729, 143 

Chen, Y. Q., Zhao, G., Nissen, P. E., Bai, G. S., & Qiu, H. M. 2003, ApJ, 591, 925 
Chomiuk, L., et al. 2011, ApJ, 743, 114 

Cresci, G., Mannucci, F., Delia Valle, M., & Maiolino, R. 2007, A&A, 462, 927 
Delahaye, F., & Pinsonneault, M. H. 2006, ApJ, 649, 529 
Denicolo, G., Terlevich, R., & Terlevich, E. 2002, MNRAS, 330, 69 
Eldridge, J. J., Izzard, R. G., & Tout, C. A. 2008, MNRAS, 384, 1109 

Epstein, C. R., Johnson, J. A., Dong, S., Udalski, A., Gould, A., & Becker, G. 2010, ApJ, 709, 447 
Fulbright, J. P., McWilliam, A., & Rich, R. M. 2007, ApJ, 661, 1152 

Georgy, C, Ekstrom, S., Meynet, G., Massey, P., Levesque, E. M., Hirschi, R., Eggenberger, P., & Maeder, A. 2012, [axXlv7l26 3.5243 
Georgy, C, Meynet, G., Walder, R., Folini, D., & Maeder, A. 2009, A&A, 502, 611 
Gunn, J. E., et al. 2006, AJ, 131, 2332 

Han, X. H., Hammer, F., Liang, Y. C, Flores, H., Rodrigues, M., Hou, J. L., & Wei, J. Y. 2010, A&A, 514, A24 
Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532 

Israelian, G., Garcia Lopez, R. J., & Rebolo, R. 1998, ApJ, 507, 805 
Jensen, A. G., & Snow, T. P. 2007, ApJ, 669, 378 

Johnson, J. A., Gal- Yam, A., Leonard, D. C, Simon, J. D., Udalski, A., & Gould, A. 2007, ApJ, 655, L33 



14 



Stoll et al. 



9.5 



9 



o 

E- 



8.5 



8 

8 8.5 9 9.5 

PP04N2 

Figure 15. Forward and reverse conversi o n ove r the entire valid conversion ranges between the met allicity scale of [ Trcmonti ct al. ( 120041) 
and the N2 diagnostic of Pettini & Pagcl (2004 ). The scales are the same as those in Figure 3 in Kcwlcv & Ellison (2008). The black 
solid line is the conve rsion defined by Kcwlcv & Ellison (2008) from t he N2 diagnostic o f|Pcttini & Pagcl (2004) to the me tallicity scale of 
Trcmonti ct al. ( 2004). The red dotted line is the conversion defined by Kcwlcv & Ellison ( 2008) from the metallicity scale of Trcmonti ct al. 
120041) to the N2 diagnostic of lPet tini & Pagcl (2004). The blue dashed line marks where a 1:1 correspondence between the two diagnostics 
would fall. The green long-dashed line shows our fit to the forward conversion, which we use in this paper to convert from T04 to PP04N2 
in order to ensure monotonic behavior at low metallicities. 

Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205 

Kauffmann, G., et al. 2003, MNRAS, 341, 33 

Kelly, P. L., & Kirshner, R. P. 2011. larXiv:1110.1377l 
Kessler, R., et al. 2009, ApJS, 185, 32 
Kewley, L. J., & Dopita, M. A. 2002, ApJS, 142, 35 
Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183 
Kobulnicky, H. A., & Kewley, L. J. 2004, ApJ, 617, 240 
Kozlowski, S., et al. 2010, ApJ, 722, 1624 

Kriek, M., van Dokkum, P. G., Labbe, I., Franx, M., Illingworth, G. D., Marchesini, D., & Quadri, R. F. 2009, ApJ, 700, 221 
Kudritzki, R., & Puis, J. 2000, ARA&A, 38, 613 

Langer, NT., Norman, C. A., de Koter, A., Vink, J. S., Cantiello, M., & Yoon, S. 2007, A&A, 475, L19 
Lara-Lopez, M. A., et al. 2010, A&A, 521, L53+ 

Lecureur, A., Hill, V., Zoccali, M., Barbuy, B., Gomez, A., Minniti, D., Ortolani, S., & Renzini, A. 2007, A&A, 465, 799 
Levesque, E. M., Kewley, L. J., Berger, E., & Jabran Zahid, H. 2010, AJ, 140, 1557 
Li, W., et al. 2011, MNRAS, 412, 1441 

Maiolino, R., Vanzi, L., Mannucci, F., Cresci, G., Ghinassi, F., & Delia Valle, M. 2002, A&A, 389, 84 

Mannucci, F., Cresci, G., Maiolino, R., Marconi, A., & Gnerucci, A. 2010, MNRAS, 408, 2115 

Mannucci, F., Salvaterra, R., & Campisi, M. A. 2011, MNRAS, 414, 1263 

Mannucci, F., et al. 2003, A&A, 401, 519 

Martini, P., et al. 2011, PASP, 123, 187 

McGaugh, S. S. 1991, ApJ, 380, 140 

McWilliam, A. 1997, ARA&A, 35, 503 

Modjaz, M., Filippenko, A. V., Silverman, J. M., Kleiser, I. K. W., & Morton, A. J. L. 2010, The Astronomer's Telegram, 2503, 1 
Modjaz, M., Kewley, L., Bloom, J. S., Filippenko, A. V., Perley, D., & Silverman, J. M. 2011, ApJ, 731, L4+ 
Modjaz, M., et al. 2008, A J, 135, 1136 

Murphy, J. W., Jennings, Z. G., Williams, B., Dalcanton, J. J., & Dolphin, A. E. 2011, ApJ, 742, L4 
Neill, J. D., et al. 2011, ApJ, 727, 15 
Niino, Y. 2011, MNRAS, 417, 567 

Ober, W. W., El Eid, M. F., & Fricke, K. J. 1983, A&A, 119, 61 

Okada, Y., Onaka, T., Miyata, T., Okamoto, Y. K., Sakon, I., Shibai, H., & Takahashi, H. 2008, ApJ, 682, 416 
Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 1979, MNRAS, 189, 95 
Pauldrach, A., Puis, J., & Kudritzki, R. P. 1986, A&A, 164, 86 
Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Prantzos, N., & Boissier, S. 2003, A&A, 406, 259 

Prieto, J., & Filippenko, A. V. 2010, Central Bureau Electronic Telegrams, 2224, 3 
Prieto, J. L., Stanek, K. Z., & Beacom, J. F. 2008, ApJ, 673, 999 
Rau, A., et al. 2009, PASP, 121, 1334 

Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329 

Reddy, B. E., Tomkin, J., Lambert, D. L., & Allende Prieto, C. 2003, MNRAS, 340, 304 

Rich, R. M., & Origlia, L. 2005, ApJ, 634, 1293 

Rich, R. M., Origlia, L., & Valenti, E. 2007, ApJ, 665, L119 

Rodriguez, M. 2002, A&A, 389, 556 




Type II SNe Environment Metallicities 



15 



Table 1 

Observation properties 



Telescope 


Instrument 


slit width 


ruling 


A coverage 


resolution 






(arcsec) 


(lines/mm) 


(A) 


(A) 


APO 


DIS 


1.5 


B400/R300 (gratings) 


3500-9800 


7 


du Pont 


WFCCD 


1.7 


400 (blue grism) 


3700-9200 


7 


Hiltner 


OSMOS 


1.2 


704 (grism) 


3960-6870 


3 


SDSS 


fiber spectrograph 


3 (dia) 


B640/R440 (grisms) 


3800-9200 


3 



Salim, S., et al. 2007, ApJS , 173, 267 

Sanders, N. E., et al. 2011. larXiy-lTTO .2363 

Smartt, S. J. 2009, ARA&A, 47, 63 

Stanek, K. Z., et al. 2006, ??jnlActa Astron., 56, 333 

Stoll, R., Prieto, J. L., Stanek, K. Z., Pogge, R. W., Szczygiel, D. M., Pojmanski, G., Antognini, J., & Yan, H. 2011, ApJ, 730, 34 
Stoll, R., et al. 2010, in Proc. SPIE, Vol. 7735, 154 
Tinsley, B. M. 1979, ApJ, 229, 1046 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 

Uomoto, A., et al. 1999, in Bulletin of the American Astronomical Society, Vol. 31, American Astronomical Society Meeting Abstracts, 
1501 

Vergani, S. D., et al. 2011, A&A, 535, A127 
Vink, J. S., & de Koter, A. 2005, A&A, 442, 587 
Wright, E. L. 2006, PASP, 118, 1711 

Yates, R. M., Kauffmann, G., & Guo, Q. 2012, MNRAS, 422, 215 

Yin, S. Y., Liang, Y. C, Hammer, F., Brinchmann, J., Zhang, B., Deng, L. C, & Flores, H. 2007, A&A, 462, 535 
York, D. G., et al. 2000, AJ, 120, 1579 
Young, D. R., et al. 2010, A&A, 512, A70 

Zaritsky, D., Kennicutt, Jr., R. C, & Huchra, J. P. 1994, ApJ, 420, 87 



Stoll et al. 



Table 2 

Measured host metallicitics 



SN Name 


z 
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12+log(0/H) 


[Fe/H] 


SN radius a ' b 
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a SN radius from the center of the galaxy. For the 12 targets for which we use archival SDSS spectra 
(indicated), this is the distance betwecnthe SN location and the fiber center, rather than from the 
center of the galaxy; a difference of lialf an arcsccond or less in all cases. 

PTF09bcl and PTF09fma are outside the SDSS photometry, so we do not determine the coordinates 
of their host galaxy centers. 

c Projected physical radius calculated (Wright 2006) assuming Hq — 70, fl-rn — 0.3, and Oa — 0.7. 
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Table 4 

Host properties measured from SDSS photometry 
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a Assuming H (i = 70, tt m = 0.3, and O a = 0.7 

Magnitudes arc corrected for galactic extinction but not intrinsic extinction. 



Type II SNe Environment Metallicities 

Table 5 

Host galaxy properties fit from SDSS photometry 
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Table 5 — Continued 
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